Universal properties of a single polymer chain in slit: Scaling versus MD simulations 
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We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic 
cispects. We confirm a number of well known scaling predictions and analyse their range of validity 
by means of comprehensive Molecular Dynamics simulations using a coarse-grained bead-spring 
model of a flexible polymer chain. 

The normal and parallel components of the average end-to-end distance, mean radius of gyration 
and their distributions, the density proflle, the force exerted on the slit walls, and the local bond 
orientation characteristics are obtained in slits of width D = 4 10 (in units of the bead radius) 
and for chain lengths A'' = 50 300. We demonstrate that a wide range of static chain properties 
in normal direction can be described quantitatively by analytic model - independent expressions in 
perfect agreement with computer experiment. In particular, the observed proflle of conflnement- 
induced bond orientation, is shown to closely match theory predictions. 

The anisotropy of conflnement is found to be manifested most dramatically in the dynamic be- 
havior of the polymer chain. We examine the relation between characteristic times for translational 
diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and nor- 
mal relaxation times are in good agreement with our observations. A novel feature is the observed 
coupling of normal and lateral modes with two vastly different relaxation times. 

We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain 
length. 

PACS numbers; 36.20-r, 36.20.Ey, 02.70.Lq 

I. INTRODUCTION 

One of the most impressive successes in the theory of polymer solutions was the discovery of a close analogy with 
critical phenomena in ferromagnetic systemsii^i^i^i^: scaling theory based on this analogy predicts that in the case of 
a good solvent and flexible polymers the variation with molecular mass and concentration of many properties directly 
measured experimentally can be presented in universal form as functions of a dimensionless parameter - the ratio of 
a characteristic length to the average size of a polymer coil (Flory radius^) Rp = aN'^ with the Flory index v k, 2>/b 
in dilute solution. Here a is the size of an effective monomeric unit and N - the number of such units in a polymer 
chain whereby, for simplicity, a prefactor of an order of unity is suppressed. For example, the mean-squared distance 
between monomers, mean-squared radius of gyration, static structure factors, were presented in scaling form^. The 
results were also generalized to include the effects of varying spatial dimensionality, solvent quality, and chain stiffness. 
Most of these predictions were verified both experimentally and by computer simulations. In fact, for an athermal 
isolated flexible macromolecule in three or two dimensions the level of description is so detailed that it includes not 
only the average characteristics but the probability distribution of the gyration radius, R, the two-point correlation 
functions and the average number of intramolecular contacts as a function of R. 
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The scaling theory was then extended to treat the effects of interfaces and geometrical restrictions. This direction 
of research is strongly motivated by important applications to colloid stabilization and flocculation, liquid chromatog- 
raphy, osmotic-pressure chromatography, ultrafiltration and others. A number of problems are thereby closely related 
to the adsorption of polymers at surfaces and interfaces. However, even in the absence of adsorption, geometric 
constraints imposed on a polymer induce dramatic changes in its behavior which have consequences pertaining to 
important applications. The basic model that captures almost all the essential physics is a single polymer chain in 
a slit or a tube with neutral (repulsive) wallsM^^ai^a3^a5a6i7^i9^,2i^,23^^^^^^^ Lately, the main 

focus has shifted to semidilute and concentrated solution o^^'^°i'^^ . Nevertheless, even the basic model has intriguing 
properties some of which have not been fully clarified yet. 

The presence of a solid impenetrable wall gives rise to depletion effects that were successfully treated within 
the scaling framework^Si^i^i^ and found experimental verification. The most dramatic effects are observed when 
long polymer chains are confined in slit-like or tube-like nanochannels where the chains become effectively two- or 
one-dimensional. This situation is typical for biological objects in living matter. Scaling theory formulated by de 
Gennes^i^ii^ describes the basic polymer characteristics, namely the average size of a chain and its confinement free 
energy. It is well understood that the scaling expressions are valid only asymptotically: R ^ D ^ a, otherwise 
corrections to scaling come into play. The scaling predictions have been successfully tested^^'^* by means of Monte- 
Carlo simulation methods^. 

In this paper we present the results of detailed Molecular Dynamics (MD) simulations using a coarse-grained bead- 
spring model with the aim of establishing the range of validity of the scaling theory for uncharged flexible polymers 
in good solvent strongly confined in a slit. One of the goals of the paper is to give simple expressions that provide 
numerically accurate values for all the essential characteristics of the chain with regard to the polymerization index, 
N, and the slit width D. 

We analyse the principle conformational parameters (end-to-end distance r and the gyration radius R) for both 
lateral and normal directions. We examine the probability distributions W{r) and W{R) for the lateral components 
of r and R, establish their universal character and provide simple analytical expressions for them. Special attention is 
payed to the properties in normal direction which have been largely ignored in earlier studies. Thus we derive analytic 
expressions for the normal components of (r^) and (i?^) which turn to be model- independent, in perfect agreement 
with the simulation data. The bond orientation profile across the slit as well as the mean bond orientation against 
slit width D have been predicted analytically and verified here by computer experiment. The average force produced 
by the polymer on the confining walls is calculated directly, suggesting thus a scaling form for the free energy of 
confinement with accurate numeric coefficients. 

An essential part of our investigation examines dynamical properties of polymers in confinement. We calculate the 
time autocorrelation functions for the gyration radius and the end-to-end distance, extract the relaxation times, and 
demonstrate distinct dynamic scaling behavior for lateral and normal relaxation aa well as their mutual interplay. 

Recently discovered possibilities of manipulating individual chains by AFM, optical tweezers and surface force appa- 
ratus have generated strong interest in the properties of end-grafted chains subject to external forces and confinement 
effects^i^. In view of these developments we present a theoretical analysis and MD simulation data on the effects of 
grafting with respect to the static and dynamic characteristics of the polymer. 

The paper is organized as follows: In Section [TT] we give a brief description of the model, and in Section IIIII we 
discuss the properties of a confined chain in equilibrium. We give a derivation of the bond orientation profile across 
the slit in the Appendix. Section IIVI examines the relaxational and diffusive dynamics of the polymer. Finally we 
discuss the implications of our results in Section |Vl 

II. MODEL AND METHOD 

The setup of our simulation is displayed in Figure [1] The polymer chains are described by a simple coarse-grained 
bead-spring model, originally proposed by Kremer and Gresli^, which has been widely and very successfully used for 
MD simulations of polymers in various contexts^i^. Effective monomers along the chain are bound together with a 
finitely-extensible non-linear elastic (FENE) potential, 

UpENEir) = -15e^(7^o/a)2 In (l - r' /TZl) , 7^o = 1.5a (1) 

where a is the range parameter of a purely repulsive Lennard- Jones (LJ) potential, that is truncated and shifted to 
zero in its minimum and acts between any pairs of monomers. 

ULj{r) = 4e [(a/r)i2 - (a/r)^ + l] , r < r, - 2^/^ (2) 

The parameter e, characterizing the strength of this potential, is chosen unity, also the temperature fc^T = 1, thus 
this potential gives rise to an excluded volume interaction betwen all the non-bonded monomers of the chain; for 
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FIG. 1: A snapshot of a tethered chain with A'^ = 200 monomers in a sht of width D — 4.0. The soft repulsion of the wall 
potential is indicated by a gradient in the color intensity. The effective width D + 5 is indicated by dashed lines. The grafting 
monomer is indicated by darker color. 



the bonded monomers, together with Ufene{'^), it leads to a typical neighbor - neighbor distance of Vmin = O.OGcr. 
Henceforth we choose cr = 1 as our unit of length. 

Similarly, the confining walls are presented also by a purely repulsive wall defined by Eq. [2] whereby the wall position 
is placed at the potential minimum. 

Molecular Dynamics (MD) simulations were performed using the standard Velocity- Verlet algorithm**^, carrying 
out typically 1.5 x 10^ time steps with an integration time step St = O.Olto where the MD time unit (t. u.) io = 
{a'^ra/ A^ei^jY^"^ = 1/V48, choosing the monomer mass m = 1. The temperature was held constant by means of a 
standard Langevin thermostat with a friction constant Co = 0-5 



III. GLOBAL EQUILIBRIUM CHARACTERISTICS 

A. Free Energy and Force 

According to the blob picture^, a chain confined in a narrow slit of width D will form a two-dimensional self-avoiding 
walk consisting of n;, blobs of size D. Each blob contains g = [pja)^!^ monomers so that the number of blobs is 
rif, = Njg = N{D /a)~^/'^ where a is the distance between neighboring monomers, and V = 0.58758(7)2^ jg Ylory 
index in 3-dimcnsional space. The free energy excess of the confined chain (in units of ksT) is simply the number nf, 
of blobs, 

F,,^, = Bnb = BN{D/a)-^/\ (3) 

where _B is a model-dependent dimensionless numerical coefficient. Note that this free energy is extensive in N . 

The confinement free energy and the associated force acting on the walls were studie d^^'^^ by means of the pruned- 
enriched Rosenbluth method (PERM algorithm) that allows one to estimate the partition function^^. For the off-lattice 
model of Ref.^-, the pressure tensor was calculated using the virial theorem and studied in conjunction with the normal 
density profile. For relatively narrow slits the simulation data for the force were consistent with the scaling prediction 

For self-avoiding walks on a simple cubic lattice the confinement free energy was calculated by MC methods in3^i^. 
These results were presented by a general power-law fit as i^conf = 0.843A^"-^'*(L)/a)~^'^^. According to this fit, the 
free energy is not exactly extensive and the I?-exponent deviates from the scaling prediction which makes the fitting 
formula purely phenomenological. Hsu and Grassberger^^ have obtained accurate data for very long chains [N up 
to 8000 and slit width up to L* = 120) and verified the scaling prediction. They also noted that there must be a 
non-universal correction to the value of the slit width D to D + 5 with 5 — 0.33, since the original scaling formula 
assumes the asymptotic limit D ^ a. Their best fit for the confinement free energy was i^conf = 2.10A^(D/a)~^/''. 
The importance of proper corrections to the slit width in order to achieve good scaling was emphasized by Milchev 
and Binder^^ and by TeraokEfc^S. 

The MD method does not allow to obtain the free energy, but the force acting on the walls can be calculated 
directly. The force acting on the walls is directly derived from the confinement free energy 



fa = -dF,oni/dD = iB/iy)N{D/a)-^-^^' 



(4) 
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Figure [2] presents the force against the scahng combination N{D/a)^^^^/'^ with and without a shift 6. It is clear 
from Figure [2] that the data points do not faU on a universal curve, if 5 = 0. The best universal fit for the force was 
obtained for S ^ a: 

fa = 2mv-^ND-^-^/''. (5) 

The corresponding free energy is given by 

Fconf = 2.03A^p/a (6) 

This correction to D turned out to be relevant not only for the free energy but for describing the variation of all other 
quantities as functions of D too. Throughout the paper we use the same slit width correction 5 = a. It is remarkable 
that the scaling asymptotic relation for the confinement free energy is achieved very early, when the number of blobs 
is still close to unity: the minimal number of blobs in our simulations was rtf, — N{D/a + l)-!/" = 0.84 for iV = 50 
and D = 10. 

Comparing our results for the free energy with those obtained for a lattice model we notice that a very good semi- 
quantitate estimate of the slit confinement free energy is 2kBT per blob. It is also quite remarkable that the same 
estimate^^ works very well for the free energy of a 2D chain confined in a strip (i^conf — 1.944n;, where nt = ND~^^'^^). 



B. Density profiles 

The distribution of the end monomer position across the slit was studied extensively by Hsu and Grassberger— . 
The scaling variable is taken as the ratio ^ = z/D and a simple scaling formula that takes into account the mirror 
symmetry with respect to the median plane was proposed. If the coordinate z is counted from the median plane, the 
expression has the following form 

The distribution Eq. [7] is properly normalized to unity. It is well known that for an ideal chain in a slit the total 
monomer density profile is proportional to the square of the end segment distribution, as long as the strong confinement 
condition is satisfied. A hypothesis that the same relation holds also for self-avoiding chains was checked by Hsu and 
Grassberger-^ and good agreement was observed. With this ansatz, the monomer density profile can be written as 



Density distributions of end-monomers and of all monomers across the slit width (for tethered chains of different 
length N and slits of different width D) are presented in Figure [3] together with the theoretical curves. The data 
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FIG. 3: Density distribution of end- monomers (a), and all monomers (b), across the slit width for tethered chains of different 
length A'' and slits of different width D (both given as parameters). Dashed lines denote the theoretical results, see Eqs.[71[S] 



points for various values of N and D collapse onto the universal theoretical curves containing no fitting parameters. 
Some deviations violating the mirror symmetry can be seen on the full density distribution picture. These are due to 
the effects of grafting which are especially noticeable for the shortest chain N = 100 in a wide slit D = 10a. The fact 
that the analytical formulas based on the scaling ansatz work very well both for the end monomer distribution and 
the full density profile for lattice as well as for off-lattice models is quite remarkable. 

The scaling ansatz has been also verified by Hsu and Grassberger— for two-dimensional chains (with 1/2 = 3/4) 
confined in a strip, and for confined ideal chains (with i^g — 1/2) in a slit. We would confidently speculate that the 
formulae will work for real flexible polymers in good and ^-solvents including exotic situations when strong adsorption 
enforces essentially 2D conformations while additional repulsive barriers confine the chain inside a strip. 

C. Average chain size characteristics. 




N(D+5) "' A^^2(^+5) wvyv) 

FIG. 4: (a) Normal components of the end-to-end vector, r±, and the gyration radius, R±, against the number of blobs 
N{D + 5)^^''". (b) Variation of the parallel components of the end-to-end vector, ry, and gyration radius, _R||, with scaling 
variable N^^iD + S)'^''"^^" for chains of length iV and different slit width 4 < D < 10. 



The characteristics for a polymer chain average size can be reduced to moments of the joint monomer-monomer 
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distribution p{zi, zj). In particular, the mean square of the gyration radius {R\) is calculated as: 



(1/27V) E / /( 



ZjYp{zi, Zj)dzidzj 



(9) 



whereby the integration in Eq. [H]is carried out across the slit. It follows naturally from the scaling blob picture of 
a chain strongly confined in a slit that the joint distribution splits into a product of single-monomer distributions 
p{zi,Zj) = p{zj)p{zi) as long as the monomers are separated along the chain by at least one blob, \j — i\ D^/'^ . We 
assume that the dominant contribution comes from the internal monomer pairs (that is, pairs which do not belong 
to the two terminal blobs). Hence the single-monomer distribution p{zj) is independent of the index j and coincides 
with the normalized full density profile p{z). It follows that 



(z — z') p(z)p{z')dzdz' 



r(2 + f) / r(3 + i) r(2 + i) 
r(i + i) ^r(4 + f) r(3 + f) 



(10) 



i. e., (-R5_) — 0.04_D^, where the scaling form of the full density was used. The derivation suggests that effects of 
grafting should not be discernable in The average square of the normal end-to-end distance is defined as 



(zi — zn)'^p{zi, ZM)dzidzM 



(11) 



We employ the same factorization assumption to separate the probability densities for the two chain ends. In the 
case of a non-grafted chain both ends are described by Potid{z) distributions leading to 



{Zi - Zn)"^ Pend{zi)pend{zN)dzidzN 



1 , 2r(2 + i) fr{3 + ^) r(2- 



2i/ ' 



2 r(i + 2^) Vr(4 + i) r(3-fi) 



(12) 



If one end is fixed (tethered) at some point Zgr, then the distribution for this end is given by p{zi) — S{zi — Zgr) 
reducing the formula to (r^(zgr)) = {r±{0)) + zL where 



^^^^^ i r(4 + i)r(i + i) 4."' 



(13) 



The result is quite general and covers a whole class of confined systems characterized by different values of the 
Flory index v. In our model the chain end was fixed at Zgr = 0/2—1 and the simulation data are to be compared 
to the theoretical expression 



. = Q.iD^ -D + l 



(14) 



The formula suggests that the combination 



D-l 



(15) 



remains constant independent of D and N . Strictly speaking, the above theoretical predictions are valid only asymp- 
totically for rif, ^ 1. The simulation data for the three characteristics calculated above are presented in Figure [Da, 
as functions of the blob number. The end-to-end distance for the grafted chains is corrected according to Eq. [151 In 
full agreement with the asymptotic theory, the normalized values are independent of D and N , and the actual values 
coincide with the theoretical predictions. Small deviations are visible only for the gyration radius curve with number 
of blobs close to 1. It is quite remarkable that the limiting asymptotic values are achieved even for such moderate 
compression. 
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It is well known that the size of a free chain in the bulk is proportional to . For the Kremer-Grest modeP^ that 
we use in our MD simulations, the model-dependent numerical prefactors defining the rms average of the gyration 
radius and the end-to-end distance have been estimated as 

Rf = 0.2707^^ and rp = 0.6707^^. (16) 

According to the blob picture, a chain in a slit is a two-dimensional self-avoiding walk consisting of nj, blobs of size 
D. Thus scaling predicts the average lateral size to be 

^nl^D^ aN"' {D/a)^-"^/" . (17) 

If D is of the same order as Rp, Eq. (fT7| . one finds R\\ ~ aN" , providing the expected smooth crossover. 

Figure [4)d presents the average lateral r|| component of the end-to-end distance and the gyration radius versus the 
scaling arguments as suggested by the theory. For the lateral size, no effect of grafting was observed. 

The best fit for the average lateral end-to-end distance is 

r|| = imaN^'iD/a + if-'^^l'' (18) 
and for the average lateral component of the gyration radius 



= 0.4007^^ (£>/« + 1)1""^/^ (19) 

Although the amplitudes are non-universal, it would have been instructive to compare them to those of other models 
in order to see the variation range. Unfortunately, only one amplitude for the end-to-end distance of chains on a cubic 
lattice r\\ = ^.%?>haN''''-{D j a + xf-'^'^l^ could be found in the literature^^. 



D. Probability distributions of end-to-end distance and radius of gyration 

The scaling form for the end-to-end vector distribution of a self-avoiding polymer chain was suggested by Fisher— 
in 1966 and then refined by Cloizauxi^ and others. In the MD simulations we obtain the histogram for the scalar 
end-to-end distance. The distribution for the reduced scalar distance r\ = {r/ro) in d dimensions can be represented 
as 



(20) 



where rp is the average end-to-end distance, and A and B are numerical coefficients. The exponent Sd = 1/(1 — i>d) 
describes strong stretching of the chain and is related to the Flory exponent Vd for a chain in the d-dimensional 
space. The exponent 6d = (7^ — l)/iyd is related to another critical exponent 7^ that appears in the partition function 
Qn = /x~^iV^''~^ where fi is called an effective coordination number. In d = 2 one has 72 = 43/33i, in d = 3 the best 
estimate is 73 ~ 1.1575(6)*-. The extra factor rQ^r]^'^~^^ that appears in Eq. [20]in comparison to the cited references 
is due to the volume element included in the definition of the distribution of the scalar distance. According to this 
definition, the two normalization conditions have the following form: J Wd{r)dr — 1 and ^ rWd(y)dr = tq, and they 
fix uniquely the vales of A and B: 
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d+l+( 
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(21) 



A 



Sd 



= 1.374 

The final form of Wd{rf) for d = 2 is 



Sd 



/r 



d + Od) 



Sd 



(22) 



1^2(77) = 1.374?7i ''58ea;p (-0.3247?") 



(23) 
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For d = 3 the distribution is given by 



Wsirj) = 3.032ry2-2^^exp(-1.08V'i25) 



(24) 



It is clear that the distribution W{ri) for a finite chain in a sht has to be between W2 and W3. For narrow shts and 
long chains, rp ^ D, the distribution W{r) is expected to be close to For wide slits or short chains with tq « D, 
W{r]) should be closer to the three-dimensional distribution W3. 
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FIG. 5: Probability distribution functions of the end-to-end distance of a polymer chain parallel to the slit wall for varying chain 
length A'' and slit width D. Lines indicate the analytic results for a confined chain (full line), Eq. 1231 and for a chain without 
geometric constraints, Eq. 1241 (dashed line). The curves are normalized so that the average value {rf) — 1 and W{{Tj)) = 1. 



In Figure [5] we plot the distribution W{ri) of the lateral component of the end-to-end distance (in semi-log scale) 
against rj = v/tq for various chain lengths N and slit widths D. The analytical equations W2 and W3 for chains 
in two-dimensions and three-dimensions are shown by solid and dotted lines, respectively. In the region around the 
maximum both distributions are close to each other and all data points collapse on a universal curve. For large 
extensions r > l.Srg, the strongly confined chains follow the two-dimensional curve W2 while for chains in relatively 
wide slits the data points lie closer to the three-dimensional curve W3. In this region W{ri) is not universal and 
depends on the ratio ro/D. The distribution of the gyration radius for a chain in a narrow tube was postulated by 
Victor^'^ in scaling form, and verified numerically by Sotta et al.^"'^ and Bishop et alliSi,. They also calculated the 
probability distribution for the end-to-end distance and fitted their results by a similar expression. The distribution of 
the gyration radius and of the end-to-end distance for a chain in a slit has not been analyzed yet. 

In Figure [5] we plot the distribution of the lateral component of the gyration radius (in semi-log scale) vs 
C, — R\^/Rq for various chain lengths N and slit widths D. The analytical equations W2 and 

W^2(C) =0.65(r^-t-C^-2) (25) 

and 

w^ic) = 1.34 (r^ + - 2) (26) 

for two-dimensional and three-dimensional chains are shown by solid and dotted lines, respectively. In the region 
around the maximum both distributions are close to each other and all data points collapse on a universal curve. 
For large extensions R > I.SRq, the strongly confined chains follow the two-dimensional curve W2(C) while for chains 
in relatively wide slits the data points lie closer to the three-dimensional curve W^3(C)- In this region W{C) is not 
universal and depends on the ratio Rq/D. 



E. Confinement effect on segment orientation 



A commonly used measure characterizing the orientation of bonds is the average value of the second Legendre 
polynomial of the azimuthal angle, (P2) = 2~^{3{coa^ 9) — 1). Experimentally this parameter appears in NMR and 
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FIG. 6: Probability distribution functions of the gyration radius parallel to the slit wall. Here curves are normalized so 
that the maximum value (max = 1 and W{C,max) ~ 1- 




FIG. 7: (a) Profile of mean bond orientation {P2{z)) = 2 ^ [3(cos^ 6} — l] across the slit where the orientation of the bonds is 
measured with respect to the Z-axis, for a real and Gaussian tethered chain. A horizontal dashed line denotes the theoretical 
prediction, {P2{z)}[([ = —j^{a/D)^, while a full line marks the scaling conjecture {P2{z)) oc (-^ — z^)^^^~^^'^\ (b) Variation 
of (P2), averaged over 2, against the effective slit width D + S. Circles denote simulation results for a Gaussian chain, squares 
- for a chain with excluded volume interactions. A dotted straight line indicates the theoretically predicted slope of —1.4, and 
a dashed line shows the theoretical result for (P2)jjj, Eg. 1431 



optical birefringence measurements. In a free non-confined coil the bond orientation is completely isotropic, and 
(P2) = 0. For a chain in a slit, it is natural to expect a preferential bond orientation along the lateral plane which 
would lead to a non-zero negative value of (-Pz). In Figure [7^ the profile of the average orientation across the slit 
obtained in the MD simulation is displayed together with the similar profile calculated for the same model but with 
the excluded volume interactions between non-neighboring monomers switched off. The values of (P2) averaged over 
all monomers irrespective of their position are plotted against the effective slit width in Figure [7|3 in a log-log scale, 
the best fit for the slopes being —2.0 for the ideal chain, and —1.44 for the 'real' chain (with excluded volume). It is 
clear that the excluded volume interactions affect the magnitude of the orientation very strongly. 

A naive estimate of the orientation effect as a function of the slit width can be obtained as follows: assume that 
the wall induces some lateral orientation only locally, within a distance of the order of a monomer size a. The local 
orientation is by itself independent of D. The fraction of monomers within distance a from the wall is given by 
p{D/2 + a) - D-^-^/" which would give D for a chain with excluded volume and D ^ for an ideal chain. Both 
estimates turn out to be well below the magnitude of the observed effect. Next, we introduce the effect of orientation 
correlations along the backbone of the chain. It is known that the simplest bond-bond orientation correlation function 
Pi for an ideal chain decreases exponentially with the distance s along the chain, exp(— ks) where k cx 1/a for flexible 
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chains. In a real chain the correlations decay according to a power law s~'^ where ui = 2(1 — v)^. We will assume 
that the correlations that transmit the orienting effect of the wall, as described by the (P2) parameter, propagate 
along the chain according to the same laws. In the spirit of the scaling theory, the contour distance s can be related 
to the normal distance from the wall, z, as s = z^/". Thus the relationship {P2{z))) oc z~^l^ can be represented in a 

normalized form as {P2{z)) oc z^)^"^/'^. Figure [7^ demonstrates that the simulation data appear in very good 

agreement with this theoretical prediction. 

By integrating over the slit width one obtains the average orientation parameter for the real chain as 

{P2) = {P2)v,aii / z-'^'"p(z)dz = const.D-^'" (27) 

J-D/2 

A similar calculation for the ideal chain gives 

{P2) = {P2)wau / e-^'p,diz)dz ^ const.D-\ (28) 

J ~D/2 

It is clear that the scaling law obtained for the real chain, (P2) ^ £)-2(i/i/-i) _ £)-i.4o agrees very well with 
the simulation data while the theoretical result for the ideal chain (P2) ^ still underestimates the observed 

effect. This suggests that apart from the local orienting effect of the wall there must be another mechanism inducing 
preferential orientation in the lateral plane. A theory describing this mechanism is given in the Appendix. The main 
result for the ideal chain is 

(P2)id = (29) 
which agrees well with the simulation data without any fitting parameters. 

IV. DYNAMIC SCALING 

There are only few simulation studies on the dynamics of a polymer chain confined in a slit. For such a chain 
only the center-of-mass diffusion has been studied by MC for an off-lattice model^. The characteristic diffusion time 
was found to scale as r ~ N^^^ which conforms to the dynamic scaling prediction r cx /N using the static scaling 
prediction R ~ N'^'^ . Here we present a study of relaxation dynamics for the gyration radius, derived from the analysis 
of the appropriate autocorrelation functions, and compare it to the center-of-mass diffusion characteristic time. The 
time-dependent autocorrelation functions were calculated as follows: 

{RjQfRitf) {R{Of? 
^^^('^^ ~ {Rm {R{QW ^ ^ 

separately for the lateral and perpendicular components of the gyration radius. 

Typical curves are presented in Figure [5) Part (a) demonstrates the effect of grafting one chain end. Evidently, this 
slows down considerably the lateral relaxation. The grafting effect on the normal relaxation is more complicated, and 
this is directly related to the shape of the ACF. The major initial portion of normal ACF is characterized by a rapid 
decay with a relaxation time unaffected by grafting. However, there is clearly a visible tail of rather small amplitude 
that is described by a much slower relaxation. A comparison with the lateral ACFs suggests that we encounter 
a classic example of weakly coupled degrees of freedom with a large difference in their time scales. The observed 
slow tails of the transverse ACFs repeat exactly the corresponding tails for the longitudinal relaxation. The semilog 
scale of the figure allows one to see the dominant relaxation time of the longitudinal ACFs although the ACFs are 
not perfectly straight and indicate some contribution from the faster modes in their initial decay. The fundamental 
longitudinal relaxation times can be estimated from the slopes of the linear part of the ACF rather reliably. The 
situation with the trasverse relaxation, however, is complicated by the fact that the initial relaxation contains also 
small-scale contributions from different modes whereas the long-term relaxation is dominated by coupling to lateral 
modes. 

Part (b) of the Figure displays the normal ACFs for various values of N and two values of Z) in a log-log scale. It 
is clear that the major initial portion of the normal relaxation is A'^-independent although it depends strongly on the 
slit width. On the other hand, the long tails demonstrate a dependence on both N and D due to coupling, as one 
would expect for the longitudinal relaxation. 
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FIG. 8: (a) Comparison of the time autocorrelation functions (ACF) of the normal and parallel gyration radius components for 
grafted and free chains with — 100 in a slit with D = 10. The intersection of C±{t) with the horizontal dotted line at 
has been used to determine the characteristic time of the fast chain relaxation perpendicular to the slit wall, shown in Fig. 8b. 
(b) Log-log plot of the ACF of the gyration radius normal component for two slit widths _D = 4 and _D = 10 and different chain 
lengths. 



The scaling theory for the longitudinal relaxation time is rather straightforward and the main ideas were discussed 
before by Descas et al41. The relaxation time scales as 



R\fr (31) 



where Cfr is the total chain friction coefficient. For a free-draining chain, Qr = -^Co with Co being the friction per 
one monomer. It follows immediately that the relaxation time is proportional to a scaling parameter 

^ ^ ^1+2^.2 + g)2{l-^.2/^^) (32) 

All the longitudinal relaxation times evaluated from the slopes of the gyration radius ACFs are plotted in Figure [S] 
against the scaling variable suggested by Eq. [32] The best fit for the relaxation times of non-grafted chains is 

r|| = 0.0053V'. (33) 

It is natural to expect that the effect of grafting one of the chain ends is dynamically equivalent to doubling the 
chain length as far as longitudinal relaxation is concerned. In this picture, the grafted end is similar to the middle 
monomer of the doubled chain whose motion relative to the center of mass contributes very little to the lateral size 
relaxation. The data points for grafted chains are described very accurately by an equation based on this picture, 



= 2^/20 0053^. (34) 



The global translational diffusion coefficient of the confined chain was obtained from the mean square displacement 
of the center of mass as a function of time. All the curves are nearly ideal straight lines. According to the Einstein- 
Smoluchovsky equation for a 2-dimensional diffusion 

((re™(0 - r,„,(0))2) = 4i?difr t = -^t (35) 

The friction coefficient in our MD simulation comes from the thermostat coupling and is set to be equal to 0.5. Indeed, 
all the simulation results are in perfect numerical agreement with the theory. It is a commonly accepted convention 
to define characteristic diffusion time as a time required to reach the mean-square displacement equal to R^. From 
this convention we obtain the following scaling fit for our model: 

Td^ff =0.02i! (36) 



Comparing Eqs.[551and[5Sl we find that the ratio of the lateral diffusion time to the time of lateral relaxation amounts 
to a constant factor of about 3.8. 
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FIG. 9: (a) Mean relaxation times parallel to the slit wall for free and tethered chains with different length A'^ in slits of width 
D against scaling variable = N^^^{D + S)'^^^~'^^^'^\ Since a tethered chain relaxes effectively as a free chain of doubled length, 
the vertical offset of the straight lines in log- log coordinates is 2^''^, as expected - see text, (b) Characteristic time of the fast 
chain relaxation perpendicular to the slit wall against width D + S for chains of different length A'^. The straight line indicates 
the theoretically expected slope of 2 + . 

A scaling description for the normal relaxation is based on the idea that, as far as trasversal motion is concerned, 
blobs relax independently. This would mean that both the normal end-to-end distance and the normal gyration radius 
are characterized by the same relaxation time of a single blob, 

As mentioned above, the shape of the transversal ACFs if far from being simple so that extracting a characteristic 
time requires caution. In order to evaluate the proper normal relaxation times we were subtracting the long-time tails 
due to coupling with longitudinal modes. After this subtraction, the shape of the ACFs allows for a better defined 
dominant transverse relaxation time. 

The transverse relaxation times are presented in Figure [Qb against the slit width D on a log-log scale for various 
chain lengths. The simulation data is consistent with the scaling prediction of Eq. [37] which gives t± oc D^+i/'', being 
is marked by a dashed best fit line in Fig.[5]3: tj_ = 0.0021?^+^/'^. However, the accuracy is not very high as indicated 
by the error bars. 

One can estimate numerically the typical time for normal relaxation, making use of Eq. 1331 by noting that T|| 
describes the relaxation of a single blob, provided one uses in ijj the number of monomers in a blob. This yields 

Tj_ = ombSD^+^Z". (38) 

which is proportional to the best fit line in Fig. [9)3 by a factor of w 2.7. 

V. SUMMARY AND DISCUSSION 

In the present work we present a comprehensive study of the static and dynamic properties of flexible polymer 
chains confined in a narrow slit with impenetrable repulsive walls. A combination of extensive MD simulations and 
an analytic theory provides a consistent picture of polymer behavior parallel and normal to the slit walls. 

• The confinement free energy (per blob) has been obtained rather precisely and shown to be equal very nearly 
to 2kBT. 

• Exact, model-independent expressions with no adjustable parameters are derived for the normal components of 
the mean end-to-end distance and radius of gyration, and shown to be in excellent agreement with simulation 
data. Also in the case of grafted chains, closed analytic expressions for the end-to-end distance in normal 
direction were obtained for arbitrary positions of the grafting monomer. 

• The observed size of the polymer parallel to slit walls is found to comply very well with scaling predictions in a 
broad interval of chain lengths and slit widths. The corresponding probability distribution functions are found 
analytically and confirmed by means of our computer experiments. 
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• The bond orientation profile across the sUt was predicted analytically for Gaussian and real chains, and verified 
by simulation. The predicted average orientation for different slit widths has been found to agree very well with 
simulation data. 

• The characteristic relaxation times of confined chains in directions parallel and normal to the slit walls have been 
obtained from evaluation of the respective autocorrelation functions. The effect of grafting on lateral relaxation 
time, which should be equivalent to doubling of the chain length, has been demonstrated. 

• A novel feature is the observed coupling of normal and lateral modes with vastly different relaxation times. 

• It is found that the mean diffusion times in lateral direction, albeit scaling similarly to the lateral relaxation 
time with chain length and slit width, are larger than the latter by a constant factor of about 3.8. The lateral 
relaxation time, evaluated for a single blob, is larger by the perpendicular relaxation time by a factor of 2.7. 

It is interesting to note that the possibility to derive accurate analytic expressions for the static chain properties 
perpendicular to the slit planes is due to the screening of correlations between distant blobs in the direction of 
compression. 

As mentioned before, the properties of the present system of a confined polymer are entirely determined by the 
underlying anisotropy in space. The situation is similar to that of adsorbed polymers on a plane where the thickness 
of the adsorbed layer is determined by the attraction to the surface and the chain conformations undergo deformation 
in perpendicular direction^. Therefore it is not surprising that the adsorbed chain dynamics is described by scaling 
theory in close analogy to the present treatment, as shown recently by Descas et al4i. 
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Appendix 

In order to evaluate an indirect orienting effect of the walls we take into account the fact that a bond experiences 
a torque due to the two tails attached to it. This torque is non-uniform as the partition functions of the tails depend 
on the position in the normal direction. The non-normalized weight for all configurations with the n-th monomer at 
position z and the n + 1-st monomer at position z' ~ z + a cos 9 is given by the product of the partition functions of 
the two tails 

W{z,9)^Q„{z)QN-n{z + acos0) (39) 

Here we assume factorization which is exact for ideal chains. (This is also a plausible scaling ansatz for real chains in 
the strong confinement limit provided only normal directions are concerned). The partition function Q{z) is in fact the 
Green's function integrated over the positions of the other tail end. Up to a normalization factor it coincides with the 
probability distribution of free chain ends pend{z), and is therefore a function of the form Q{z) = constD~^ f {z / D) . 
For a confined ideal chain in the ground state approximation it is independent of n and is given by the ground state 
wavefunction of the Edwards' equation: 

Q(z) = const. D~'^ cos{ttz/D) (40) 
Expanding Q(z + a cos 6) up to second order in cos 9 we obtain 



W{z,9)^Q'{z) 



l + -jcos0+-(-j -cos 9 



(41) 



Calculating the average (cos^ 6*) with the weight W up to the lowest non-trivial order in a/D gives the following 
result for the average Legendre polynomial 



(42) 
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In the general case this quantity is still a function of the z coordinate and has to be averaged over the slit width 
with the weight given by the full monomer density p{z). However, for an ideal chain, the ratio /"// = """^ is constant 
at any position within the slit and won't be affected by averaging. This brings the final result 

(P2>id = -^ (43) 

which means that the indirect orientation effect is dominant compared to the local wall effect ^ D^^ 

For a real chain, after averaging over the slit width we get a result very similar to Eq. 1431 Only the numerical 
coefficient changes slightly but the dependence persists. Since the propagating effect of the wall is much stronger 
in this case, the indirect orientation by torque constitutes just a minor correction. 
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